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A formalism is presented for analytically obtaining the probability density function, Pn{s), for 
the random distance s between two random points in an n-dimensional spherical object of radius 
R. Our formalism allows P„(s) to be calculated for a spherical n-ball having an arbitrary volume 
density, and reproduces the well-known results for the case of uniform density. The results find 
applications in stochastic geometry, computational science, molecular biological systems, statistical 
physics, astrophysics, condensed matter physics, nuclear physics, and elementary particle physics. 
. . . As one application of these results, we propose a new statistical method obtained from our formalism 

04 , to study random number generators in n-dimensions used in Monte Carlo simulations. 
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I. INTRODUCTION 



(N 

(N. „^ . . . 

In two recent papers [0, 0], geometric probability techniques were developed to calculate the functions P3 (s) which 
• describe the probability density of finding a random distance s separating two random points distributed in a uniform 
> ' sphere and in a uniform ellipsoid. As discussed in Refs. [0, |, |, |, |, ||, |, |, |l^, |l|, |l2|, |l4[ |l| , these results 

■ are of interest as tools in mathematical physics, and have numerous applications in other fields as well. Specifically, it 
I was demonstrated in Refs. ||l|, ^, ||, [To[ |l^ that knowing the random distance distribution in a spherical object greatly 
. facilitates the calculation of self-energies for spherical matter distributions arising from electromagnetic, gravitational, 

■ or weak interactions. 

. As an example we calculate the total electrostatic energy W3 of a collection of Z charges uniformly distributed 
' within the same spherical volume of radius R. For illustrative purposes, we assume that Z is a large number. For 
each pair of charges the potential energy due to the Coulomb interaction in Gaussian units is 

- — - (1) 



where e is the elementary charge, ri (7^2) is the coordinate of the first (second) charge. The total Coulomb energy 14^3 
can then be expressed as 

Wz = \^ j ^i^^i j sin^id^i j d(j)i j rldr2 J sin 02^6*2 J — d4>2 



s 

> 

• I— I 

C3 ■ = — X , (2) 

^ ■ 2 5 i? ' ^ ' 

where p = SZe/AirR^ @. We note that Eq. (|) requires evaluating a six-dimensional integral, and using the additional 
theorem for spherical harmonics , 

^'^'°'^) = 27TT ^ (-iri^™(^i,'^i)>^r'"(^2>2). (3) 

m— — l 

Alternatively we can use the probability density function P3 (s) giving the random distance distribution for a sphere 
with a uniform density to calculate 1^3. For a collection of Z charges there are Z [Z — 1) /2 such pairs, and hence the 
total Coulomb energy W3 is: 

W3 = — ^ X / P3 (s) — ds = — X 9^ — X , (4) 

^ 2 Jn s 2 5 i? 2 5 i? ' ^ ' 
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where ^ 

We note that by using geometric probabihty techniques we can simpUfy the expression for W3 from a six-dimensional 
non-trivial integral (Eq. (|^)) to a one-dimensional elementary integral (Eq. (Q)). Generalizing to n-dimensions, a 
calculation of the electrostatic energy Wn for a collection of Z charges uniformly distributed within the same n- 
dimensional spherical volume of radius R, can be greatly simplified by utilizing the n-dimensional random distance 
distribution. This reduces the complexity of calculating Wn from a 2ri-dimensional intractable integral involving 
n-dimensional spherical harmonics, to a simple 1-dimensional integral. 

The probability density function Pri(s) for the distribution of the random distance s between two random points 
in a uniform spherical n-ball is well known. Hence the object of the present paper is to generalize the results of 
Refs. § 0, 1 to the case of an arbitrarily non-uniform density distribution by using a new method which we present 
below. Notice that the sample space B„ for the random points is a spherical n-ball of radius R defined as 

Bn = {{xi,X2,--- ,Xn) e R'' : xl + xj + ■ ■ ■ + xl < R^} , (6) 

where 

RJ" ^ {x : X ^ {xi,X2r ■ ■ ,Xn)} (7) 

represents n-dimensional Euclidean space. 

To illustrate our formalism, we begin by deriving the probability density function (PDF) for a spherical n-ball with 
a uniform density distribution, and compare our results to those obtained earlier by other means [^,0, 1). We then 
extend this technique to a spherical n-ball with an arbitrary density distribution, and this leads to a general-purpose 
master formula for P„ (s) given in Eq. (^9|). The outline of this paper is as follows. In Sec. || we present our formalism 



and illustrate it by rederiving the well-known results for a circle and for a sphere of uniform density. In Sec. Ill 



we extend this formalism to the case of non-uniform but spherically symmetric density. In Sec. [V we develop the 
formalism for the most general case of an arbitrary density distribution. In Sec. ^ we present some applications of our 
results. These include the mth moment (s™) for a spherical space with a uniform and Gaussian density distribution, 
the evaluation of the self-energy arising from i/z/-exchange interactions, the probability density functions for a sphere 
with multiple shells that arises in neutron star model jl^, and finally a new proposed computational scheme for 
testing random number generators in n-dimcnsions. 



II. UNIFORM DENSITY DISTRIBUTIONS 



In this section we illustrate our formalism by deriving the PDF for a circle of radius R having a spatially uniform 
density characterized by a density function p, where p is an arbitrary constant. For two points randomly sampled 
inside the circle located at fi and r2 measured from the center, define a random vector s = r 2 — ^1 and a random 
distance s = \s\, where < s < 2R. To simplify the discussion, we translate the center of the circle to the origin so 
that the equation for the circle is x'^ +y'^ = R^. It is sufficient to initially consider those vectors s which are aligned 
in the positive x direction, since rotational symmetry can eventually be used to extend our results to those vectors 
s with arbitrary orientations. We begin by identifying those pairs of points, fi and r2, which satisfy s = sx. One 
set of random points for fi is located in Ai and the other set of random points for r2 is located in A2 as shown in 
Fig. |l|. We observe that A2 is the overlap area between the original circle Ci and an identical circle C2 whose center is 
shifted from (0, 0) to (|s | , 0) as shown in Fig. ||. Since the areas of Ai and A2 are equal, it follows that the probability 
density of finding a given s = |s2:| in a circle of uniform density is proportional to the area of A2: 

I dx I dy+ I dx I dy. (8) 

Js-R J-y/R:^-{x~sf Jf J-VR2-x^ 

Using rotational symmetry, this result can apply to any orientation of s, where < (j) < 2tt, and hence the probability 
density P2 (s) for a circle of uniform density can be factorized as 

P2 (5) = 27rs X / (s) , (9) 

where / (s) is a function to be determined. If we impose the normalization requirement 



p2R 

/ P2is)ds^l, 
Jo 



(10) 
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FIG. 1: Locus of points in a circle separated by a vector s — sx. For each random point ri in Ai. there is a unique random 
point r2 in A2 such that s = f2~ri. 




FIG. 2: The probabihty density that two random points are separated by a random distance \sx\ in a circle of radius R is 
proportional to the shaded area given by the overlap of Ci and C2, where Ci is x + y = R and C2 is (x — s) + y — R . 



we then have 



2s 

i?2 



dy) 



ds 



4s 



-i(s/2i?). 



(11) 



(12) 



Equation (|T^ is identical to the results obtained in Refs. ||, Q by other means. 

The conclusion that emerges from this formaUsm is that the probability density of finding two random points 
separated by a random vector s in a circle of uniform density can be derived by simply calculating the overlap region 
of that circle with an identical circle obtained by shifting the center from the origin to s. In the following discussion, 
we show that this result generalizes to higher dimensions, and provides a simple way of calculating P„(s) for n > 3. 

The above formalism for a circle of uniform density can be extended to a sphere of uniform density. For a given s, 
we select the positive z direction and study the distribution of the random vectors s in this direction. It then follows 
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that the probability density of finding sz is proportional to 

dz 



s-R 



dx 



dy - 



dz 



dx 



dy. 



(13) 



In a 3-dimensional space, we note that the direction of the random vector s can have the following range: < 9 < it 
and < 4> < 27r. Following the previous discussion, we thus arrive at the following expression for a sphere with a 
uniform density distribution: 



ins' 



{l^dzl 



-vR^-z'^ "-^ j-vR^-z^-x^ "-y 



/o'' 4^«' {l" dzj^^j^ dxj^^^^^^ dy^ ds 
„ s2 9 s3 3 s5 



(14) 
(15) 



i?3 4 i?4 16 i?6 ■ 

The result in Eq. ( p^ ) agrees exactly with the expression obtained previously in Refs. ^, ^, 0, ||]. 

The present formalism can be readily generalized to express P„ (s) for a sphere of uniform density in n-dimcnsions 

as 



Pn{s) 



s" ^ f^,„ dxn dx-i ■ ■ ■ f 

Js/2 " J_^R2_^2 i J_ 



-y/R'^-xl----xl_ 



dx 



n-l 



s«-i f" dx r 

Jo I * Js/2 "-^"J.^ 

We find that if n is an even number, 



^ - dxi ■ ■ ■ J 



'-sjR'^-x\----x\ 



(16) 



— COS 
TT 



-M^/2i?)--E 



fe=i 



(n-2fc)!! 
(n-2fc + l)!! 



V4)^ 



(17) 



where 0! = 0!! = 1. If ti is an odd number. 



Pn{s) = 



R" 



1)!! ^ 



(n-l)!! f-' 2fc + lfc! ( 



n-l 
2 



fc)! 



1 - (s/2i?) 



2fe+l 



(18) 



The functional forms of Eq. ( |l6| ) for n = 1, 2, 3, and 4 are shown in Fig. ||. 
The cumulative distribution function (CDF) |l], ^ for Pn{s) is given by 



D„(x) 







P„(s)ds 

R /'I 2 _|_ I"! 
_ -Pg V2 ' 2 2/ 

^(i t + 5) 



^ + 2" X ( 2 2 ' 2 2 ^ 



1 H 4_ 1^ 
2' 2 2; 



(19) 



where < x < 2i?, a = x^/4i?^, and i?ct is the incomplete beta function. 

We summarize three important representations for the probability density function P„ (s) for a spherical n-ball of 
radius R with a uniform density distribution as follows: 



1. Integral representation: 



Puis) 



dx 



i,i)i?2" 



(20) 



2. Generating function representation: 



Pn(s) 



sin ^ h — sin" 



-1 / 2/i-\/4^ 
V2-hV4^ 



2n-">. 2 2' 2/ 



(21) 



where \h\ < 1 and i? = 1. We note that one can obtain a number of identities and recursion relations for P„ (s) 
from the generating function representation in Eq. (|2l|). 
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3. Hypergeometric function representation: 



n ('n I 1 1^ 
^ V 2 2 ' 2^ 



2F1 (a, b,c,a)R~ - 2F1 (a, 6, c, /?) 



(22) 



where a = 1/2, 6 = 1/2 — n/2, c = 3/2, a = 1, /3 = s^/4i?^, and 2^1 (■ • • ) is the hypergeometric function [ p7[ . 
Notice that we can obtain the orthogonahty relations for P„ (s) from the hypergeometric function representation 
as shown in Eq. (p^. 



In Rcf. |Tl| additional identities and recursion relations for P„ (s) are discussed in greater detail. 



III. SPHERICALLY SYMMETRIC DENSITY DISTRIBUTIONS 



In this section we extend the previous results to the case of a circle with a variable (but spherically symmetric) 
density characterized by a density function p (r) . Following the derivation presented in the previous section, we note 
that for any random vector s = sx, if the second random point r2 carries the density information p{x,y), then the 
first random point fi should have the density information p{x — s,y) . It then follows ^2(5) can be expressed as 

s dx J^^T^ p {x ~ s, y)xp {x, v) dy 

' r(-tf-/-^^<^--)-<-)*)- 

A general formula for P„ (s) for an 71-dimensional spherical ball of radius R having a spherically symmetric density 
can be derived, and we find 



S" dXn J ^ r— dxi--- j I— 5 5 pi X p2 dXn-1 



2 4 
S 



2 4 
S 




0.5 r 
0.4 

CD 

0.2 
0.1 





FIG. 4: Plots of Ps (s), P4 (s), -P5 (s), and Pe (s) for a Gaussian density distribution. Note that in all cases the functions are 
normalized such that Pn{s) ds = 1. 



where 



Pi = pixi,X2, ■■■ ,x.n- s), 

P2 = p{xi,X2, ■ ■■ ,Xn) ■ 



(25) 



Some analytical results for a sphere with various spherically symmetric density distributions can be found in Ref. pi] . 

As an application of Eq. (|24| ) we consider the case of an n-dimensional spherical space of radius R ^ 00 with a 
Gaussian density given by 



where 



Pn{r) 



N — lim n- 



N _i .2 
e 2^ 



(27r) 2 cr 



Pn{r)r" ^dr. 



(26) 



(27) 



In Eq. (g^) r is measured from the center, and the integral is over all space. The PDF for an 71-dimensional spherical 
Gaussian space can then be obtained as 



Pn{s) 



2"-lr (f) cr"' 



(28) 



Fig. ^ displays the functions in Eq. ( |2^ ) for n = 3, 4, 5, and 6, where cr = 1. Finally, we note that the maximum 
probability density, denoted by Smax, occurs at 



yj2{n^l)a. 



(29) 
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IV. ARBITRARY DENSITY DISTRIBUTIONS 

We consider in this section the probabihty density functions for a spherical n-ball having an arbitrary density 
characterized by a density function p. We begin with a circle and use the conventional notation for polar coordinates: 
X = r cos (f> and y = r sin (f>. In a 2-dimensional space, the random vector s can be characterized by an angle in the 
range < 4> < 27r. Associate each random unit vector s {(j)) with a rotation operator ^ such that 

(30) 

To ensure that the product of p (ri) and p (r2) maintains the correct density information, we use a 2 x 2 matrix 



-R2x2(0) 



cos (p — sm ( 
sin 6 cos 6 



(31) 



to characterize this particular operator such that 

p{x, y) — > p (cos (l)x — sin (j) y, sin (j)x + cos (j> y) . (32) 

Notice that R2x2{4>) is an orthogonal matrix which satisfies R2x2i^) ~ -^2x2{'t') and its determinant is +1, where T 
denotes the transpose. 

We can then express P2 (s) for a circle with an arbitrary density distribution as 

s C 1^/2 I^^w^ P {x', y') X P ix", y") dy 
r'2(s) — — ^ r (33) 

where 



p{x',y') — yO (cos (/) (a; — s) — sin0 y, sin (a; — s) + cos 0?/) , 
p{x",y") — p (cos (f> X ~ sin y , sin 4> X -\- cos (f) y) . 



Figure || exhibits ^2(5) when R — 1, and illustrates the agreement between the Monte Carlo simulation and the 
analytical result. 

The preceding discussion can be generalized to the case of a spherical n-ball with an arbitrary density. Define the 
following hyperspherical coordinates 

xi = rsin0„_2sin0„_3 sin 02 sin 0i cos (/i, 

X2 = r sin9n-2sin9n-3 sin 6*2 sin 6'i sin 0, 

Xi = r sin6'„_2 sin0„_3 sin0i_i cos6'i_2, 

x„-_i = r sin6'„_2 cos6'„_3, 

Xn = rcos6'„_2- (34) 

where 9i € [0, tt] and 4> G [0, 2tt]. Associate each random unit vector s{9i, - ■ ■ ,0) with a rotation operator ^ such 
that 

\x„) =I%\s{9i,--- ,cj))). (35) 
The n X n matrix representation for the rotation operator ^ in Eq. ( |35| ) can be expressed as 

Rn 9i, 4>) — Rnxn{4>) X Rnxn{(^l) X • • • X i?„ 

), (36) 
where n > S. The various matrices on the right-hand side of Eq. ( |36| ) are defined as follows: 

1- Rnxni4')' The matrix elements are an = cost/), 012 = — sint/i, 021 = sincp, 022 — cos^, for i € [1, 2] and j €E [3, n] 
ttij = 0, for i G [3, n] and i € [1, 2] — 0, for i G [3, n] and G [3, n] aij — Sij. 
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2. Rnxn {Oi)- The matrix elements are an = cos^i, ai2 — 0, ai3 = sin^i, 021 = 0, 022 = 1, 023 = 0, 031 = — sin^i, 
032 = 0, 033 — cos6'i, for i G [1,3] and j € [4, n] — 0, for i € [4, n] and j G [1,3] a,y = 0, for i G [4, n] and 
j G [4,n] ay = ^ij. 

3. Rnxn (Ok) and A: G [2, n- 3]: The matrix elements are afc+i,fe+i = cos^fc, afc+i^fc+2 = sin0/c, 0^+2,^+1 = - sin^^, 
afc+2,fc+2 = cos 9k, for i G [1, A;] and j G [1, k] Uij = Sij , for i G [fc + 3, 71] and j G [fc + 3, n] Oij — Sij , for i G [1, /c] 
and J G [fc + 1, n] aij = 0, for i G [fc + 1, n] and j G [1, fc] fly = 0, for i G [fc + 1, fc + 2] and j G [fc + 3, n] aij = 0, 
for i G [fc + 3, n] and j e [fc + 1, fc + 2] = 0. 

4. Rnxni0n-2)- The matrix elements are a„_i_„_i = cos6'„_2, an-i,n = sin6'„_2, a„,n_i = -sin6'„_2, a„,„ = 
cos 9n-2, for i G [1, n — 2] and j G [1, — 2] aij = Sij , for i £ [l,n—2] and j G [n — 1, n] a^j = 0, for i G [n — 1, n] 
and J G [1, 7^ — 2] a^j = 0. 

Notice that all the matrices in Eq. ( ^ ) are orthogonal and their determinants arc +1. The matrix elements for 
Rnxn{0n-2, ' ' ' ,^1,0) Can bc Summarized as follows: 

1. The matrix elements for the 1st colmrin are an — cos 9i cos (j), 021 — cos 0i sin(^, 031 ~ — sin^i, and for i G [4, n] 
an = 0. 

2. The matrix elements for the 2nd column are ai2 = — sin0, 022 = coscf), and for i G [3, n] 0^2 — 0. 

3. The matrix elements for the jth column where j G [3,n — 1] are aj+ij = — sin0j_i, for i G [j + 2,n] = 0, 
and for i G [1, j] Ojj = cos0j_i x x^, where Xi = f x Xi and is the ith Cartesian coordinate component for a 
unit vector f in j-dimensions. For example: 

ai4 = cos 03 sin 02 sin 01 cos 0, 
024 = cos 03 sin 02 sin 0i sin 0, 
034 — cos 03 sin 02 cos 0i , 

044 = cos 03 cos 02- 

4. The matrix elements for the nth column are for i G [1, n] a.m = Xi, where Xi = f x Xi and Xi is the ith Cartesian 
coordinate component for a unit vector f in n-dimensions. 
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Additionally, it is convenient to define the following transforms: 

n 

x'i = ^aij {xj -6jns), (37) 



x'- = ^flyXj, (38) 



where are the matrix elements for Rnxn 7 • • • , ^i, defined in Eq. (|3^). The master equation of P„(s) for a 
spherical n-ball with an arbitrary density characterized by a density function p (x) = p [xi, X2^ • • ■ , Xn) can then be 
formulated as 

s""^ X F„ (x,e) , , 

Pn{s) = ^^^^ , 39 



where 



Fn{x,9) = / sin"-2 0„_2 d0«-2 / sin"-3 0„„3 
Jo Jo 



n-3 «t7n~3 ' ' ' 



sin 9i dOi I dcj) dxn / dxi 

Jo Js/2 J-y/R'^-xl 



P {x') X p {x") dXn-1, 



and 



p{x') ^ p{x\,X2,-- • ,0 , 

p{x") = pixixi--- ,x:). 

The probability density function for the random distance distribution including Euclidean distance and geodesic 
distance for an n-dimensional sphere (i.e. the boundary of a spherical n-ball) with an arbitrary surface density 
distribution will be discussed in Ref . . 

V. APPLICATIONS 

A. mth Moment of the Distance 

In some applications the mth moment of the distance, rather than the distance itself, is of interest. As an example, 
for a collection of nucleons interacting via simple harmonic oscillator potentials, (s^) may be of interest rather than 
(s) itself. We then calculate the mth moment (s™) for the case of a spherical n-ball of uniform density, where 

(s") = / 5"P„(s)ds. (40) 
Jo 



Using Eq. (EOf), the mth moment (s™) has the general form 



.n+m I \ ^ V 2 ' 2 ' 2 ' 2 ' 2 



V2~'~2'2''"2"'~ 2 ) 



■R"\ (41) 



where B (p, q) is the beta function, and m > — (n — 1). 

Additionally, (s™) can be evaluated for a spherical space having a Gaussian density distribution where p cx 
From Eq. (BS) we find. 



r2B. p ( n+m \ 

\ 2 I 



= lim / s'"P„(s)ds = (2cr) 



r(t) 



(42) 
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In some applications (such as in nuclear physics) involving low-energy interactions among nucleons the lower limit 
(zero) should be replaced by the hard-core radius = 0.5 x 10"^'^ cm In such cases the expressions for P„(s) 
and (s™) can be expressed as follows: 



Pn{s) 



and 



where 



and 



H {R, Tc, m, n) 



C (a, m, n) = 



(2i?)"^ 



2R 



s"-i/f (^'-2;2)'^ dx 
C{2R,0,n) - C{rc,0,n) 

H{R, rc, m, n) 



H{R, r„0,n) 



m+7i — 1 



ds 



/2 



n -\- m 

n+m 
_ j_c 

n + m 



„ , n In 1 m\ „ ( n In 1 m 

S - + -.-H \- — ] - {- + -,- -\ \ 

2 22 2 2/ l2 2'2 2 2 



B 



1 n 1 



2' 2 2 



Bx 



In 1 



2' 2 2 



(43) 



(44) 



(45) 



with X — r1/AR^, and m an integer. We note that Eqs. ( ^3|) and ( ^4|) are the first known analytical results for Pn{s) 
and (s™) to incorporate the hard-core radius r^. 



B. Neutrino-Pair Exchange Interactions 

A second example of interest is the i/i/-exchange (neutrino-pair exchange) contribution to the self energy of a nucleus 
or a neutron star. For two point masses the 2-body potential energy is given by H, |2^, |2^, 28, 2^, 30, 



Vuy (\ri - rj\) 



47r3| 



(46) 



where and aj are coupling constants which characterize the strength of the neutrino coupling to fcrmions i and j 
(«, j — electron, proton, or neutron). In the standard model ]32|, 



2 sin^ 0w = 0.964 



2 sin^ Ow = 0.036 



1 

As an example, we find for the case of a uniform density distribution of radius R containing N neutrons, 



N{N - 1) 
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3r. \ G 



2 \2r^R3 ArcR'^ 8R^ 16R<^ J Att^' 
where is the hard-core radius. The analogous result for a Gaussian density distribution is 



J.2 4^2 



N{N_-l)Gj 
32tT37r7/2 



(47) 



(48) 



where r(a, 5) is the incomplete gamma function. 
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C. Neutron Star Models 

Another application of current interest is the self-energy of neutron star arising from the exchange of vv pairs p9| , 
|30| , Here we evaluate the probability density functions in 3 dimensions for neutron stars with a multiple-shell 
density distribution, which is what is typically assumed in neutron star models |l^, |T^. For illustrative purposes, we 
discuss a spherically symmetric model with 2 spherical shells, each of uniform density, where for simplicity we assume 
shells of equal thickness. Some other multiple-shell models and their n-dimensional probability density functions can 
be found in ||ll|. 

For a 2-shcll model with a uniform density in each shell define p = pi for < r < R/2 and p = P2 for R/2 < r < R, 
where pi and p2 are arbitrary constants and r is measured from the center of the neutron star. Using the preceding 
formalism we can show that ^3(5) has 4 different functional forms specified by 4 regions: 



1. < s < ii? 



^ ^ 24{pl + 7pl)s' _ SGjpl - 2p,p2 + 5pl)s^ ^ 12{pl - 2p,p2 + 2pl)s^ 
ipi + 7p2fR' ipi + 7p2fR^ {pi + lp2fR^ ' 



2. \R<s<R 



81{pi-p2)p2S , 24pis2 36pi(/9i +3p2)s3 



3. i? < s < |i? : 



2{pi + 7p2fR^ (Pi + 7p2)i?3 (pi + 7p2)'i?^ 

Upjs^ 
'(pi + 7p2)'i?6' 



81(P1 -P2)/52S , 24(9pi - P2)/02S2 36(5pi - P2)P2S^ 



(50) 



2 ifl + 7p2f i?2 (pi + 7p2f i?3 (pi -I- 7p2)' 

12(2pi - P2)P2S^ 



(Pl+7p2)'i?6 



(51) 



4. f i? < s < 2i? 



Psis) = T T \ 9 . (52) 

ipi + 7p2fR^ {pi + 7p2fR^ {pi + 7p2fR^ 

We observe that the probability density functions defined in adjacent regions are continuous across the boundaries 
separating the regions. 



D. RIPS: A New Test for Random Number Generators in n-Dimensions 



Another interesting application of the present work is as a new test of random number generators in n-dimensions. 
Our statistical method follows by applying the probability density functions for the random distance distribution to 
evaluate the expectation values of fi2 ■ r23 in ri-dimensions, where fi2 = T2 — r\, r23 = ^3 ~ '^2, and ri, r2, and fa are 
three random points independently sampled from a spherical n-ball. The quantity (ri2 ■ ?^23)„ is one of the geometric 
probability constants ||ll| and has applications in many areas of physics . It follows from the preceding formalism 
that for a spherical n-ball of radius R with a uniform density distribution, 

(ri2-f23)„ = -^i?^. (53) 

Note that 



lim (ri2 • f23) = -i?^. 

n — >oo 



(54) 
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TABLE I: Comparison of random number generators (RNG) with the exact results derived in this paper for (ri2 • '^23)„- The 
number of simulation samples in each case is A'^ = 10^. We note that the NWS generator fails to pass the new randomness test 
in all dimensions that we selected. See text for further discussion. 



RNG 


n — 3 


Result 


n — 5 


Result 


n = 10 


Result 


RANO 
R31 

NWS 


-.60037 ± .00065 
-.60028 ± .00066 
-.64119 ± .00071 


Pass 
Pass 
Fail 


-.71483 ±.00059 
-.71468 ±.00058 
-.75427 ± .00057 


Pass 
Pass 
Fail 


-.83310 ± .00047 
-.83229 ± .00048 
-.85016 ± .00049 


Pass 
Pass 
Fail 


Exact 


-.60000 




-.71429 




-.83333 





2 2 

In a spherical Gaussian space where p oc e"^ ' and R —> oo, 



{ri2 ■ r23)n 



(55) 



We refer to this test as RIPS (Random Inner Product in a Sphere). 

We now apply the RIPS test to check three popular random number generators frequently used in Monte Carlo 
simulations |33, Isl, E^, Is^, Is^. The random number generators tested are: 



1. RANO IQ which is a linear congruential generator and uses the following algorithm 

/„ = 16807 X In-i mod m = 2^^ - 1. 

2. R31 ||3^ which uses the generalized feedback shift register (GFSR) method 

where p — 31, q ~ 3, and ® is the bitwise exclusive OR operator. 

3. NWS |3^, |3^ which uses the nested Weyl sequence 

Y„ = {n {na}} , 

where {y} is the fractional part of y and a is an irrational number. 



(56) 



(57) 



(58) 



The results are shown in Table |. We have checked several initialization methods and the results are not affected. 
We note that the NWS generator using the nested Weyl sequence method fails to pass our randomness tests in all 
dimensions selected. Hence caution should be exercised in using this particular random number generator for Monte 
Carlo simulations, especially in molecular dynamics simulations. More statistical tests for a variety of random number 
generators including the PSLQ algorithm using the binary digits of tt [p8| , and other newly proposed computational 
schemes for n-dimensions, will be discussed in greater detail in Ref. |39|. 

The geometric probability constants, (fi2 • ^23)„, can be added to the family of various computational tests for 
random number generators. They can then serve to investigate the quality of random number generators for questions 
of randomness, especially in higher dimensions (n > 3) where few results are currently available. The possibility that 
some random number generators such as NWS pass other tests but not ours, may indicate that the properties of 
random points in a sphere provide a more sensitive test of randomness than is otherwise available. These and other 
issues will be discussed in more detail in Ref. iKOl . 



VI. CONCLUSIONS 



A formalism has been presented in this paper for evaluating the analytical probability density function of the random 
distance distribution for a spherical n-ball with an arbitrary density distribution. We show that the random distance 
distribution technique can reduce otherwise difhcult calculations from the complexity of 2rt-dimensional integrals 
to just a 1-dimensional integral, even when the nuclcon-nucleon hard-core radius Tc is included. Our formalism has 
applications to the currently active area of research surrounding string-inspired theories of higher dimensional physics, 
and has numerous potential applications to other fields as well. Specifically the results presented here are of interest 
in the context of recent work on the modifications to the Newtonian inverse-square law arising from the existence of 
extra spatial dimensions ||40|| . We have also presented a new computational method to test random number generators 
in n-dimcnsions used in Monte Carlo simulations. 



13 



VII. ACKNOWLEDGMENTS 



The authors wish to thank Chinh Le, Michelle Parry, David Schleef, Dave Seaman, Christopher Tong, and A. W. 
Overhauser for helpful discussions and the Purdue University Computing Center for computing support. This work 
was supported in part by the U.S. Department of Energy under Contract No. DE-AC02-76ER01428. 



D. Schleef, M. Parry, S. J. Tu, B. Woodahl, and E. Fischbach, J. Math. Phys. 40, 1103 (1999). 
M. Parry and E. Fischbach, J. Math. Phys. 41, 2417 (2000). 

R. Deltheil, Ann. Fac. Sci. Toulouse 11, 1 (1919). 
J. Hammersley, Ann. Math. Stat. 21, 447 (1950). 
R. Lord, Ann. Math. Stat. 25, 794 (1954). 

M. G. Kendall and P. A. P. Moran, Geometrical Probability (Hafner, New York, 1963). 

L. A. Santalo, Integral Geometry and Geometric Probability (Addison- Wesley, Reading, MA, 1976). 

A. W. Overhauser, unpublished. 

E. Fischbach, Ann. Phys. 247, 213 (1996). 

E. Fischbach, D. E. Krause, C. Talmadge, and D. Tadic, Phys. Rev. D 52, 5417 (1995). 

S. J. Tu, Ph.D. thesis, Purdue University, West Lafayette, Indiana (2001). 

N. Fedotov, Pattern Recognition and Image Analysis 8, 264 (1998). 

K. Jansons, Proceedings of the Royal Society of London Series A 432, 495 (1991). 

K. Mecke, International Journal of Modern Physics B 12, 861 (1998). 

R. Miles, Computer Graphics and Image Processing 12, 1 (1980). 

J. D. Jackson, Classical Electrodynamics (Wiley, New York, 1999), 3rd ed. 

G. Arfken and H. Weber, Mathematical Methods for Physicists (Academic Press, San Diego, California, 1995), pp. 749-750, 
4th ed. 

D. Pines, J. de Physique 41, 111 (1980). 

S. Shapiro and S. Teukolsky, Black Holes, White Dwarfs and Neutron Stars (Wiley, New York, 1983). 
W. Feller, An Introduction to Probability Theory and Its Application (Wiley, New York, 1967), 3rd ed. 

C. W. Gardiner, Handbook of Stochastic Methods (Springer- Verlag, New York, 1985), 2nd ed. 

N. G. van Kampen, Stochastic Processes in Physics and Chemistry (North-Holland, Amsterdam, 1992). 
W. Fleming, Functions of Several Variables (Springer- Verlag, New York, 1977). 
S. J. Tu, to be published. 

A. Bohr and B. R. Mottelson, Nuclear Structure Vol. 1 (W.A. Benjamin, New York, 1969). 
G. Feinberg and J. Sucher, Phys. Rev. 166, 1638 (1968). 

G. Feinberg, J. Sucher, and C. K. Au, Phys. Rep. 180, 83 (1989). 
S. D. H. Hsu and P. Sikivie, Phys. Rev. D 49, 4951 (1994). 

B. Woodahl and E. Fischbach, Phys. Rev. D 61, 125010 (2000). 
B. Woodahl and E. Fischbach, Phys. Rev. D 62, 045025 (2000). 
B. Woodahl and E. Fischbach, Phys. Rev. D 63, 065006 (2001). 
Partical Data Group, Eur. Phys. J. C15, 1 (2000). 

D. E. Knuth, The Art of Computer Programming Volume 2 (Addison- Wesley, Reading, MA, 1981). 

W. H. Press, S. A. Teukolsky, W. T. Vetterling, and B. P. Flannery, Numerical Recipes (Cambridge University Press, New 
York, 1992). 

I. Vattulainen, T. Ala-Nissila, and K. Kanlaala, Phys. Rev. Lett. 73, 2513 (1994). 
B. L. Hohan, O. E. Percus, T. T. Warnock, and P. A. Whitlock, Phys. Rev. E 50, 1607 (1994). 
K. V. Tretiakov and K. W. Woici echowski, Phvs. Rev. E 60, 7626 (1999). 



D. Bailey and H. Ferguson, URL tittp : //www. nersc. gov/news/pi_random072401 .html 
S. J. Tu, to be published. 

N. Arkani-Hamed, S. Dimopoulos, and G. DvaU, Phys. Rev. D 59, 86004 (1999). 



